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Abstract. 

Polarised Semi-Inclusive Deep Inelastic Scattering (SIDIS) processes, £{Se)p{S) IhX 
allow to study Transverse Momentum Dependent partonic distributions (TMDs), which reveal 
a non trivial three dimensional internal structure of the hadrons in momentum space. One of the 
most representative of the TMDs is the so-called Sivers function that describes the distribution 
of unpolarized quarks inside a transversely polarized proton. We present a novel extraction 
of the Sivers distribution functions from the most recent experimental data of HERMES 
and COMPASS experiments. Using suitable parametrizations, within the TMD factorization 
scheme, and a simple fitting strategy, we also perform a preliminary exploration of the role of 
the proton sea quarks. 



The Sivers distribution function A fg/p^{x, k±) embodies the correlation between the nucleon 
spin and the quark intrinsic transverse momentum and gives the number density of unpolarized 
quarks q with intrinsic transverse momentum k± inside a transversely polarized proton p^, with 
three- momentum P and spin polarization vector S, 

fq/pt{x, = fg/p{x, /c_l) + ^ A^fg/ptix, A;_l) S • (P x fc_L) . (1) 

fq/p{x,k±) is the unpolarized x and k± dependent parton distribution, and the mixed product 

S ■ {P X k±) gives the parity invariant azimuthal dependence. The distribution of quarks 
in polarized hadrons is not axially symmetric, and shows a non trivial 3-D structure. The 
conceptual importance of the Sivers function, its relation with the color QCD dynamics and its 
gauge properties, are being more and more understood. Dedicated experiments are running, 
under construction or planned. 

In Ref. [Ij, we presented an extraction of the Sivers distribution functions based on a fit of 
SIDIS experimental data from the HERMES and COMPASS collaborations. At that time, the 
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Figure 1. The Sivers single spin asymmetry, A^j^ , as a function of x for SIDIS production 

of vr^ (left panels) and tt~ (right panels) at HERMES. The upper (lower) panels show the results 
obtained from a fit which includes only valence (valence and sea) quark contributions. 



only direct evidence of a non-zero Sivers effect came from HERMES experimental data on the 
azimuthal moment A^J^'^'^ measured for pion and kaon production off a proton target [2j; 
the COMPASS collaboration, instead, found a Sivers effect compatible with zero for SIDIS off a 
deuteron target [3]. The extraction of the Sivers functions was made even more challenging by 
the fact that the ^'J^^'^''"'^^) measured by HERMES for K+ was surprisingly large, about twice 
as much as the analogous asymmetry for vr"*" production, which hinted at a possible important 
role of the sea quarks [I] . 

Since then, new experimental results have become available: COMPASS have released new 
data on SIDIS production of spinless hadrons off a proton target, showing a clear Sivers 
asymmetry jl]; a new HERMES data analysis, based on a much larger statistics, while confirming 
the previous pion data, softens the enhanced peak in the Sivers azimuthal moment [Sj. 

It is then timely and useful to perform a new fit, trying to understand the most essential 
features of the Sivers distribution function as shown by the present (scarce) experimental data. 

The SIDIS transverse single spin asymmetry (SSA) Afj''^'^'' "^^^ measured by HERMES and 
COMPASS, in the 'j* — p cm. frame and at order k±/Q, is given by [6l[7]: 
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(ps and (j)h are the azimuthal angles identifying the directions of the proton spin S and of the 
momentum of the outgoing hadron h respectively, while ip defines the direction of the incoming 
(and outgoing) quark transverse momentum, k± = (cos sine/?, 0); "^g^ — is the unpolarized 

cross section for the elementary scattering Iq — )• Iq. Finally, D^{z,p±) is the fragmentation 
function describing the hadronization of the final quark q into the detected hadron h with 
momentum Ph', h carries, with respect to the fragmenting quark, a light-cone momentum 
fraction z and a transverse momentum pj^ . 

The Sivers function is parameterized in terms of the unpolarized distribution function, as in 



Figure 2. The Sivers single spin asymmetry, A^jjj^ ^ , as a function of x for SIDIS production 
of (left panels) and K~ (right panels) at HERMES. The upper (lower) panels show the results 
obtained from a fit which includes only valence (valence plus sea) quark contributions. 



Ref. in the following factorized form: 

A%/pt{x,k^) = 2Mq{x)h{k^)f^,p{xM) , (3) 

with 

where Nq, Ug, (3q and Mi (GeV/c) are free parameters to be determined by fitting the 
experimental data. We adopt a Gaussian factorization for the unpolarized distribution and 
fragmentation functions with the Gaussian widths {k\) and {p\) fixed to the values found 
in Ref. [7] by analysing the Cahn effect in unpolarized SIDIS: {k\) = 0.25 (GeV/c)^ and 
{p\) = 0.20 (GeV/c)^. For the unpolarized, A:_L-integrated distribution and fragmentation 
functions we use the GRV98 [8] and DSS [9] sets. 

We best fit the HERMES proton and COMPASS deuteron data from Refs. [3l[5] and exclude 
the COMPASS proton data [3] due to the presence of some experimental scale uncertainty which 
is difficult to account for in the present analysis. The new forthcoming data from COMPASS 
on a proton target will become very useful in the future. 

In order to evaluate the significance of the sea-quark Sivers contributions we first perform a 
fit of the SIDIS data including only Sivers functions for u and d quarks. The results we obtain 
are rather satisfactory, with a total xlo/ °f about 1.05. They are shown in the upper panels 
of Figs. [1] - [4| The corresponding Sivers functions are plotted in the left panel of Fig. [5] It 
is interesting to notice that the Mi parameter, which fixes the Gaussian width of the Sivers 
function (i.e. its distribution in k±_) is rather well constrained and turns out to be between one 
half and two thirds of the unpolarized distribution function width. 

We then perform a second fit, in which we allow for u, d, s, s sea contributions to the Sivers 
function: overall, we obtain similar results, shown in the lower panels of Figs. [T] - |4j but the 
total Xdof improves to about 0.9. The corresponding Sivers functions are plotted in the right 
panel of Fig. [5] 

At this stage, our results can only be considered as qualitative and preliminary, as the 
experimental data are not yet so well established and abundant to allow a precise definite 




Figure 3. The Sivers single spin asymmetry A^J^ ^ , as a function of x for SIDIS 7r+ (left 
panels) and vr" (right panels) production at COMPASS on a deuteron target. The upper (lower) 
panels show the results obtained from a fit which includes only valence (valence plus sea) quark 
contributions. 



determination of the Sivers sea distribution functions. We can only try, by comparing the upper 
and lower panels of Figs. [T]-|4j to draw some general conclusions. 

Adding the sea quark contributions to the Sivers functions does not induce any variation in the 
^sin(0h ci>s) g^^imuthal moments for vr"*" and production, which are confirmed to be essentially 
generated by the fragmentation of a valence u quark from the proton; instead, the production of 
7r~ and mostly K~ are more sensitive to the presence of the sea quark contributions, and might 
offer the chance to gain some preliminary information on the polarized proton sea. 

We also learn that the A^J^'^'^ '^^^ asymmetries, as delivered by the presently available 
experimental data, are not very sensitive to the strange quark Sivers distributions, while they 
feel the presence of u and d contributions more strongly: we have actually checked that the Xdof 
is not altered when adding or removing s and s. 

When comparing the left and right plots in Fig. [Sj one notices that, while the u, s and s 
Sivers functions are consistent with zero, the d Sivers function can be rather large and negative, 
definitely not consistent with zero. This originates from the small-x behavior of the COMPASS 
data for tt~^ production off a deuteron target, most sensitive to d: as these data show a negative 
trend at small x (and very tiny error bars. Fig. |3]), the fit naturally tends to assign a negative 
value to the d Sivers function. 

Further experimental measurements from present and future facilities are strongly needed to 
reach a better knowledge on the important Sivers distribution functions, in particular in the 
large and small x regions where they are presently largely unconstrained. Some help might also 
come from the existing large SSAs in pp ^ hX processes, which can be interpreted within a 
TMD factorization scheme [lOl ITT] . 
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Figure 4. The Sivers single spin asymmetry ^ , as a function of x for SIDIS 

(left panels) and (right panels) production at COIVIPASS on a deuteron target. The upper 
(lower) panels show the results obtained from a fit which includes only valence (valence plus 
sea) quark contributions. 




Figure 5. First A;_L-moment of the Sivers distribution function, xA f^^^^, as a function of x 

as obtained by fitting SIDIS experimental data from HERMES and COIVIPASS. The left (right) 
panel shows the results obtained from a fit which includes only valence (valence plus sea) quark 
contributions. 



